Stability of Relativistic Matter With Magnetic Fields 
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Stability of matter with Coulomb forces has been proved 
for non-relativistic dynamics, including arbitrarily large mag- 
netic fields, and for relativistic dynamics without magnetic 
fields. In both cases stability requires that the fine structure 
constant a be not too large. It was unclear what would hap- 
pen for both relativistic dynamics and magnetic fields, or even 
how to formulate the problem clearly. We show that the use 
of the Dirac operator allows both effects, provided the filled 
negative energy 'sea' is defined properly. The use of the free 
Dirac operator to define the negative levels leads to catastro- 
phe for any a, but the use of the Dirac operator with magnetic 
field leads to stability. 

PACS numbers: 03.65.-w, ll.lO.-z, 12.20.-m 



Since the work of Dyson and Lenard in 1967 |]1| it has 
been well understood that matter consisting of N elec- 
trons and K static (but arbitrarily positioned) , point nu- 
clei, and governed by nonrelativistic quantum mechan- 
ics, is stable. (See also and see [|| for reviews.) 
This means that the ground state energy is finite and is 
bounded below by a universal constant times the number 
of particles, i.e., £^0 > const. (N + K). If the nuclei are 
either non-static or non-point the problem becomes eas- 
ier. The reason the nuclei are taken to be static is that 
they are so massive compared to the electron that if the 
nuclear mass (or the nuclear radius of 10^ ^'^cm, for that 
matter) played any important role then ordinary mat- 
ter would have to look very different from what it does; 
electron orbital radii would be many orders of magnitude 
smaller than the Bohr radius. 

We know that atoms, molecules and bulk matter exist 
and have ground states that are well described by the 
Schrodinger equation with infinitely massive, point nu- 
clei. This description needs corrections, however, even at 
the low energies appropriate to atomic structure. Effects 
due to special relativity and interactions of electrons with 
(quantized or unquantized) magnetic fields are signifi- 
cant. Usually, one appeals to quantum electrodynamics 
(QED) for these corrections, but that theory, as usually 
formulated, is not known to exist outside of perturbation 
theory. More importantly, the effects mentioned are often 
not perturbative, despite the fact that the fine structure 
constant a = /Tic ~ 1/137 is small. It is well known, 
for instance that the Dirac equation breaks down when 
Za [Z = nuclear charge), the quantity of greatest im- 
portance, is 1. The quantity Za is not at all small for 
Uranium and hence relativistic and electrodynamic ef- 



fects are not small. It is sometimes said that some proper 
high energy theory (e.g., QCD or strings) will straighten 
things out, but it is difficult to believe that we cannot 
have a theory of chemistry without theoretical insertions 
in the GeV range. 

With this in mind, attempts have been made over the 
years to approach a theory of matter by ever increas- 
ing generalizations of the many-body Schrodinger equa- 
tion, all the while maintaining the many-body aspect and 
a non-perturbative, rigorous treatment. One significant 
step was the introduction of a 'relativistic' kinetic en- 
ergy into the Schrodinger equation by replacing /2ra 
by \/ + m? ~ rac? (with p = —ihSI as usual) or, more 
simply, byp = |p|. We can also include a (non-quantized) 
magnetic field B(x) = curlA(x) acting on the orbital 
motion, i.e., we use |p + A(x)|. The Hamiltonian (in 
appropriate units, in which all dimensional parameters 
reside in a and p = — iV) is then 



N 



iJ,.ei = ^|p, + A(x,)|+aV; 



(1) 



where the Coulomb potential of the electrons (with co- 
ordinates Xi) and nuclei (with coordinates Rj) is 

JV K 

i=l j=l \<.i<j<N 

+ ^ |R,;-R,ri. 

(It is purely for convenience that we take the nuclei to 
have the same Z.) 

For A EE 0, Conlon g showed that matter was stable 
in this theory if a is 'small enough'. The constants were 
improved in |^ and finally in Q stability was proved all 
the way up to the critical value appropriate to a 'hydro- 
genic atom' (which is Za < 2/7r in this theory j|] 
instead of Za < 1 as in the Dirac theory) if and only if 
a itself is less than some critical value ac (independent 
of Z). In 0] it was shown that ac > 1/94. 

For arbitrary, nonzero A stability was shown to hold 
for (nl) in [no[ if a < a, where 



1/5= (-)Z + 2.80 



1.30 



(2) 



which permits Z < 59 for a = 1/137. Eq. (|) will be 
useful later. 
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Another important ingredient was the inclusion of ar- 
bitrarily large magnetic fields and, at the same time the 
replacement of p^/2m by the non-relativistic Pauli oper- 
ator (in suitable units) 

Ta = • (p + A(x))]2 = (p + A(x))2 + ^ . B. (3) 

This field, A(x), can be viewed in two ways: either 
as an applied external field or as some sort of average 
(in the sense of Hartrce theory) of the field produced by 
all the electrons themselves. The replacement of by 
[(pH- A(x)]^ is a generalization that has long been known 
to cause no problems with stability in the non-relativistic 
case. The essential new ingredient in the Pauli operator 
is the Zeeman term cr • B, which results in the fact that 
there is no lower bound to the energy, even for hydrogen 
JjH/ . (We have in mind here that we make no restriction 
on the field A, which is allowed to be arbitrarily large.) 
To remedy the instability due to arbitrarily large fields 
we add the field energy 

i/fieid = [8na]-^ J B{^f d^x 

in the units we have been employing. 

For a single atom it was shown a decade ago that the 
addition of -fffloid stabilizes the system {Eq is finite) if 
and only if Za^ is less than some comfortably large, but 
finite value The 'only if part relies on a deep result 

Stability with arbitrary N and K and kinetic energy 
(^ was first solved by Fefferman (unpublished) for suf- 
ficiently small a and Z; see for an announcement. 
Shortly thereafter a simple proof with good constants 
(e.g., we can take Z — 1050 when a = 1/137) was given 
p5[ but, as in the 'relativistic case' mentioned above 
there is an etc such that the system is unstable when- 
ever a > ac, no matter how small Z may be. Using the 
developments in ||l^ , it was shown in that this non- 
relativistic theory remains stable when the magnetic field 
is quantized, provided an ultraviolet cutoff is introduced. 

Thus, two essential modifications of the Schrodinger 
energy (replacing by p or the inclusion of spin- 
magnetic field interaction) require a bound on a, in addi- 
tion to a bound on Za or Za^, for many-body stability. 
One might surmise that the combination of the two would 
lead to further difficulties — even, possibly, instability for 
all a. This is not so, as we show in this letter. Indeed, 
our result is a bit surprising for it reveals an unantici- 
pated delicacy about the manner of defining an electron 
in Dirac's theory. 

The best way to include relativity and spin-field inter- 
action is to use the Dirac operator 

Va^\^ + 4(x) + I3m. 

As usual, = a ■ O, where and f3 are the 4x4 
Dirac matrices. Our many-body Hamiltonian (with no 



pair production) is then, formally, 

N 

-f^Dirac = ^ T^Aii) + CtVc + -fffiold- 

Unlike the usual (or the Pauli) kinetic energy opera- 
tors, the Dirac operator is not bounded below, so we use 
Dirac's prescription of filling the negative energy 'sea'. 
Concretely, this means that all the electrons lie in the 
positive spectral subspace of the Dirac operator, i.e., if 
A"*" denotes the projector onto the positive spectral sub- 
space for all the electrons then we allow only many-body 
functions rp (which are antisymmetric functions of space 
and four component spinors) that satisfy 

A+V = V'- 

Another way to say this is that the ground state energy 
is given by the infimum of the variational quantity (with 
respect to tp and A{x)) 

£{^P) = (V'lA+FDiracA+l^) /(A+^IA+V-). 

The idea of studying this problem goes back to Brown 
and Ravenhall who studied a single Helium atom 
with A = and A_|_ being the positive spectral subspace 
for the free Dirac operator, Vq. This problem with a 
single general atom, i.e., K = 1, but N and Z arbitrary, 
has recently been shown stable, if Za < 2/{tt/2 + 2/tt), 
and instable otherwise [ p^ . 

While the Hamiltonian H does not depend on the def- 
inition of the positive spectral subspace, and it always 
includes the A field, the meaning of A+, and hence Eq, 
does depend on the subspace. 

There's the rub! Which Dirac operator should we use 
to define A+? There are at least two significant choices. 
One is the free Dirac operator Vq = -I- (3m, which is 
the usual choice. The other is the Dirac operator with 
the magnetic field Pa • The latter seems more important 
since an electron can never get rid of its own magnetic 
field. Moreover, the choice Vq, is not gauge invariant, 
i.e., the multiplication of an electron wave function by 
a spatially varying phase usually takes a positive energy 
function into a mixture of positive and negative energy 
functions. The second choice is manifestly gauge invari- 
ant. 

This issue is usually not clearly stated in field the- 
ory textbooks, but whatever one might think about the 
appropriateness of either definition the interesting and 
surprising fact is that "stability" of matter can settle the 
argument. The following two theorems are the main re- 
sults in ||l9| . They show that one choice is valid and the 
other is not. 

THEOREM 1 (Stability with the magnetic Dirac 
operator): 
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If the electron wave function is restricted to lie in the 
positive spectral subspace of the Dirac operator Va then 
^^Dirac is stttblc (i.e., -ffDirac > 0) providcd a and Z are 
small enough. In particular, if we define ac to be the 
unique solution to the equation (with a in (Qj and Li 3 
defined below) 

1 - (ac/a)' = (le^Li^gae)"/' 

then a < ac suffices for stability. In particular, Z ~ 56 
and a = 1/137, or Z — 1 and a — 1/9 are sufficient for 
stability. 

We shall outline the proof of this theorem here; details 
can be found in ||l9t] . 

THEOREM 2 (Instability with the free Dirac op- 
erator): // the electron wave function is restricted to lie 

in the positive spectral subspace of the free Dirac operator 
Vq then the Hamiltonian -ffDirac is unstable. More pre- 
cisely, for any given values of a > and Z > 0, there are 
(sufficiently large) particle numbers N and K for which 
the infimum of £{ijj) is ~oo. 

Remark on the proof of Theorem 2: This requires a 
complicated construction of variational functions with ar- 
bitrarily negative energy. It is necessary to construct iV- 
particle antisymmetric trial functions, ijj, of space-spin 
and also corresponding magnetic potentials A(x). See 

Proof of Theorem 1: Using and (j^) we have that 



Vr>- 



1 ^ 

^ V|p, + A(x,)| . 

rv ^ — ^ 



4=1 



and hence f ("0) — -ffficid is bounded below by the sum of 
the negative eigenvalues of the one-body operator 

h := A+f2?A -k|p + A(x)|)a+, 

-ff field ^ — Tr/i_, where h_ = 

~: ifx > 0, r > 



with K — a/a. I.e., £{'4') 
[\h\-h]/2. 

Next we note the BKS inequality 
are self-adjoint operators, then 



Tr[X - y]- < Tr[X^ - Y^]_! 



211/2 



(4) 



We use (|) with X = A+VaA+ and Y = A+k\p + 
A(x)|A+. But we note that 2?^ Ta. + > Ta. > (p + 
A)2 - B(x) and we note that X"^ = K+V\A+ (since Va 
commutes with A+), while Y'^ < K^A+(p -I- A(x))^A+. 
Thus we conclude that for every choice of A and ip 



£{■4) - iJfieid > -Tr (1 - k2)(p + A(x))2 - B(x) 



1/2 



by the Lieb-Thirring inequality for the sum of the square 
roots of the negative eigenvalues of a Schrodinger opera- 
tor (p -I- A(x))^ + V(x.) with arbitrary A and arbitrary 
potential V. This is 



Tr 



(p + A(x))2 + l/(x) 



1/2 



< Li 



y^(x) d^x . (6) 



Note that in (^) there is a factor of 4 on the right side 
because the trace includes a trace over 4-dimensional 
spmors; m (|) there is no such trace, i.e., (||) is for spin- 
less particles. It is known that Li ^^ < 0.06003 . 

The factor 4 in (|) can be reduced to 2, as one might 
expect physically; details are in Theorem 1 then 

follows easily from (||). 
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> -4(1 - n^)-^/^Lx^3 / B(x)2 d^x (5) 
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